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ABSTRACT 

We use the canonical quantization of spherically symmetric dust universes, 
and calculate the expectation values of the quantized metric, < >. 
Though the classical solutions are singular, and the wave functions have no 
zero support on singular geometries, the expectation value (for a general 
quantum state) is everywhere regular. For a quasi-classical (coherent) state, 
| a >, the expectation value, < a\g\a > describes a universe (or star) that 
collapses to a minimum radius, the Planck radius, and re-expands again. 
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The canonical quantization of simple dust universes is known for many 
years [1], and it was shown recently that the quantum states describing those 
systems are bound states [2,3]. In this work we discuss the issue of singular- 
ities of those models. 

We use the geometrical coordinates, G — c — 1, and H ^ 1. 

Consider first a homogeneous close universe. It can be described by the 
Friedmann metric with the line element 



ds 2 = -N 2 ± (t)dt 2 + a 2 (t)dn 2 3 



(1) 



where dfl^ is the volume element on the 3-sphere. The matter in this universe 
is dust, with the Lagrangian [3,4] L D = —8n J y/^gpU^U u d 3 x. Where p is 
the dust density. The total Lagrangian including the gravitational part, can 
be written as [2,3] 
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where a$ is the maximum radius of the universe. The Hamiltonian of this 
system is 



H = N± 
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where P a = dL/da. The classical solution of this simple system is the well- 
known cycloidal solution [5] 



a (t) = ^°o(l + cos(f)) 



(4) 



where we chose the gauge in which N(t) = a(t). The solution (4) describes a 
universe, starting with a maximum radius, a , that collapse to the singularity, 
a = 0. 

The quantization of this system is elementary, using the Hamiltonian (3), 
we can get the Wheeler-DeWitt equation [2,3] (or the Schrodinger equation) 
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where we have used the coordinate representation < a|\I/ >= ip(a) , a = 
a , P a = —ihd/da. If we define x = a — |a , we get the harmonic oscillator 
equation 

n 2 <p 1 



+ -moj 2 x 2 
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where m = 12tt 2 , uj = 1 and E = 3n 2 al. We see that the quantum state 
describing this universe is a bound state, ip n {x), with the spectrum 

E n = 37r 2 (a 2 )(n) = h(n + ^j (7) 

We see that the initial radius of the universe can take only discrete values 
[2], 

aoW = -^V /n+1 /2 (8) 

The Oppenheimer-Snyder (OS) collapsing star is a spherically symmetric 
homogeneous dust star that collapses to form a Schwarzschild black hole [6]. 
The interior is a slice of the homogeneous Friedmann universe. One can use 
(1), but the angle x is limited to the region [0, xo\- Where Xo is a constant 
angle smaller than ir/2. The star's surface is at x — Xo- Outside the star we 
have the (static) Schwarzschild space-time. The matching conditions on the 
star's surface give 

2M = aosin 3 xo and R = a sinxo (9) 

where R is the initial Schwarzschild radius of the collapsing star, and M is 
its mass. 

The Lagrangian in that case is L = L inside + L outside . Where L inside 
is (2) (where the x~ integration is from zero to xo om y)- The outside La- 
grangian describes an empty Schwarzschild space-time, and when we study 
the semiclassical states of this system, it does not contribute. This is be- 
cause ip WKB ~ exp(iS c i ass /h), and S , c z a ss(Sch.) = 0. The classical solution in 
that case is again (4). It describes a star, initially with a radius Rq, which 
collapses to the Schwarzschild singularity, a = 0. 
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The Wheeler-DeWitt equation in this case, is very similar to (5), with 
the only difference that now < x < Xo- The states are bound states, and 
the quantization conditions are [3] 



where ao = \2ti J xo sin 2 dx- From (9) and (10) we see that the mass and 
initial radius of the Oppenheimer-Snyder star are subject to a quantization 



One way of avoiding the classical singularity is if the wave function has 
zero support on singular geometries. Namely, if VKflV = 9av ng ^) = 0- I n 
our case the singular geometries are the ones for which a = 0. This is the 
"big crash" in the cycloidal solution, and the black hole singularity in the OS 
model. If we had ip(a = 0) = 0, then the probability of having a singularity 
is zero, V S i ngu i ar = \tp(a = 0)| 2 = 0. Unfortunately, this is not the case in our 
model. The wave functions that satisfy the Wheeler-DeWitt equation are 
the standard harmonic oscillator wave functions, and those do not vanish for 
a = 0, or x = —\clq- Does this mean that we have singularities also in the 
quantum model? In other words, is the condition ip(a = 0) = 0, necessary 
for avoiding the classical singularity? 

We will suggest a much weaker condition, based on the assumption that 
the basic physical quantity is the metric, g^ u . In that case, in the quantum 
theory, if the system is in a state |\& >, one should calculate the metric 
expectation value < ^1^^]^ >. This is a real matrix. The properties of 
this real matrix determine the geometric structure of the "quantum system" , 
especially, the singular structure of it. Using this weaker condition, we will 
see that the above models can be "quantum mechanically singularity free" . 

Excplicitly, we propose the following criteria for a "quantum singularity" : 
Given a quantum state of the system, >, one should calculate the expecta- 
tion value of the metric < ^\g^ u \^ >. The quantum system is singular only 
if its expectation value is. 

The metric can be singular for some states but not for others. So this 
criteria may serve as a "wave function selector". The stronger criteria, ip(a = 
0) = 0, is not satisfied by any state in the Hilbert space of the solutions of (6). 




(10) 



[3]. 
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We are going to see that the weaker one, regularity of < g^ v >, is satisfied 
by a general quantum state. 

Consider first the "most classical situation". This is described by the 
quasi-classical coherent state, \a >, which satisfies (see the solution (4)) 

d x \a{t) >= a(t)\a(t) >= —^=a e~ tuJt \a > (11) 

where d x = \{X + P x ). Now let us calculate < a\g^ u \a >. Using (1) we see 
that the metric (in the gauge N(t) = a(t)) is 

ds 2 = -a 2 (t)dt 2 + a 2 (t)dx 2 + a 2 ! (t)sin 2 ''xdttj (12) 

and we need to calculate only < a|a 2 |a >. Where a is the operator X — \a§. 
This is elementary 

< a\a 2 \a >=-(<* + a*) 2 + -^L (a + a*) + — + - (13) 
ii 2 v i ^ v y 42 v/ 

and using (11), we get 

a 2 h 
< a\d 2 \a >= ^(l + cos(t)) 2 + - (14) 

The first term on the r.h.s of (14) is exactly the classical solution (4), and 
the second term is the quantum correction. The quantum correction depends 
of course on the quantization scheme that we use. Namely, it depends on 
the ordering that we choose. To get (13) we use the "symmetric ordering", 
namely, X 2 = d x d x + d x d x + a 2 , + (aj.) 2 . If we define the p-ordering such that 
X 2 = (1 — p)d x d x + pd x d x + a 2 + (aj.) 2 , we see that the quantum correction 
goes like ph. So for example if we use the Wick order, p = 0, we get zero 
quantum correction. 

From (1) we see that a singular quantum metric is a one in which the 
expectation value < \I/ 1 a, 2 1 \I/ >, vanishes for some t. But according to (14) 
this is impossible^. The expectation value < a|(?|a: > is everywhere regular, it 
does not vanish for any t. Remember that \a > is the semiclassical (coherent) 
state, and so one may suspect that the quantum corrections are at least as 

2 Unless p = 0. 



5 



strong for a general quantum state. This is indeed the case. The Hilbert space 
of these models is spanned by {\n >}, the harmonic oscillator eigenstates. 
For a general eigenstate, \n>, we have (using the p-order) 

< n\X 2 \n >= hn + ph (15) 

the first term on the r.h.s of (15), hn, is the one that contributes to the 
"classical value" of X, while the second term is the quantum correction. We 
see that it is the same for any \n>. Any state |\& > is a (normalized) linear 
combination of \n >-states (like \at >). So for any state |\& >, we will get the 
same quantum correction. 

The expectation value (14) describes a semiclassical dust universe, start- 
ing with a maximum radius tig, and collapsing to a minimum radius^ 

drain ~ k = lp (16) 

here ~ is because we should have the p-factor in (14). This semiclassical 
universe does not reach the classical singularity, a = 0. The dust density, at 

^ Q"min IS 

Pmax ~ -3— = Nip 3 (17) 
^min 

where N is the number of dust particles ("baryons"). 

According to our assumption that the metric is the basic physical quan- 
tity, all other physical quantities should be derived from gj%) =< a\g uu \a >. 



,(«) 

So we consider gffi as if it is a classical metric. For example the corresponding 
Ricci scalar is 



R[g [a)] = — 



ag(l + cos(t)) 2 + 21% 



-2 
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So the "quantum Ricci scalar" is bounded from above, R m ax ~ (a 

In the case of the OS model, the solution (14) describes a dust star that 
collapses to a minimum radius, lp, and does not reach the Schwarzschild 
singularity. 

These results can be easily generalized to the (spherically symmetric) in- 
homogeneous case. Spherically symmetric dust matter in 4D Einstein gravity 

3 In the geometrical units Mp — lp = h 1 ^ 2 . 
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is an integrable system [7]. One can fix the (diffeomorphism) gauge com- 
pletely, and get an infinite set of independent (non propagating) degrees of 
freedom [1]. So we will first fix the gauge classically, and then quantized the 
reduced system. One can fix the gauge by choosing the comooving coordi- 
nates in which the metric has the form 



where p and A are functions of r and TZ, and dQ 2 . is the volume element in 
S 2 . Explicitly, we choose r and TZ such that 



where U v are the 4-velocities of the dust particles, and prime denote TZ- 
differentiation. This gauge fixing corresponds to iVj_ = 1 , N l = , and it is 
complete [1]. One then gets the reduced Hamiltonian [1,3] 



where p = yl — TZ 2 , and p s is the surface of the dust star. Y(r,p) is the 
(reduced) field variable, Y = 8e A and P Y = 5L/5Y = 3ttYY /{64p^l - p 2 ). 
A most important feature of (19) is that there are no p-derivatives in it. This 
means that the (infinite number of) degrees of freedom are decoupled. One 
can descretize (19). Define 



ds l = _ dr 2 + e 2 W 2 + e 2X dn 2 



(18) 




(19) 



k 



k = 1,2, ...,7V 



(20) 



Pk 




and write (19) as a sum 




(21) 



where Y k (r) = Y(r,p k ), P k (r) = P(r,p k ), and 
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The classical equations of motion are Y = 5H re d/ 5Py and Py = —5H re d/SY. 
Given the Cauchy data, Y(p,r = 0) , P Y (p,r = 0) (or {F fc (0) , P fc (0) }), we get 
a unique solution to the classical equations [7]. 

The Hamiltonian (19) (or (21)) is time independent so we can write it as 

N 

Hred = ^Ek (22) 
k=l 

where 

E ^^{ a ^ + ^) < 23 > 

and yf = n(r = 0), Pf = P„(t = 0). 

To quantize this system, we use the coordinates representation 

Y = Y ; P Y = -ih-^ (24) 

and get the Schrodinger equation (the Wheeler-DeWitt equation) 

where we used the U YP Y ordering" (the field Y is always to the left of its 
conjugate momenta). One can use the discrete form of (25) 

Using the fact that we have a decoupling, we can write the quantum state, 
as a direct product 

|* >= |^i > |*2 > • • -\^n > (27) 

and from (22), (25) and (26) we get a set of N independent (Schrodinger) 
equations 

7f i~ ak Y k ^ + PkYk ) = Ek ^ k{Yk) (28) 



which is the equation for an harmonic oscillator [3], with mass, m k = l/2a&, 
frequency, u k = Sp k , and energy, e k = El/4N 2 j3 k . So the spectrum of (26) is 
quantized, and we can write >= \n>k >• The Hilbert space of solutions 
of the Schrodinger equation (26) is spanned by 

{\V >} = {\m > \n 2 > ■ ■ -\n N >} (29) 

The quantization conditions, e k = hu k (n k + 1/2), give 

2 _ 1 12nhp 2 k 
h k - T7 / A n k + 1/2) 



and using (22) we get 



H„d = ^ 1/2 ^ E M"h + V2) 1/2 (30) 
k 

where A k = 4v / 37rpfc/(l — p|) 1//4 . For a finite N, the spectrum (30) is discrete. 

In the case of a dust star that collapses to form a Schwarzschild black 
hole, H red is proportional to the mass of the star. In that case (30) is the 
"generalized Bekenstein area quantization" [8,3]. 

Because A k is of order 1, an average n k is of order H^ ed /h, which for 
astronomical objects is a huge number^. So we are in the semiclassical regime. 

Consider a semiclassical collapse, starting at rest (at r = 0). In that case 
we have Pjf^ = 0, and (28) is now 



h d 



2 



' ak Y k W! + (3k{Yk - Yk 



DefiningX fc = n-^ (0) , we see that {X k , P k } is a one dimensional harmonic 

oscillator. A quasi-classical (coherent) state, corresponding to a classical 

(cl) 



(Oh 



<f k (Y k ) = (31) 



solution X k , is 



x (o) 

K >= -^e-^> fc > (32) 



V2 



1 For our sun for example, we have Mq/% ~ 10 6 
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We want to calculate the metric expectation values. These are < \l/|e 2Al |\I/ > 

and < \I/|e 2A |\[ / > (see (18)). In our gauge however, we have e M = ^ 1 %p P dY/dp. 
So we need to calculate only < ak\Y£\ak >■ Again, this is elementary, we 
use the p-ordering and get 

< a k \Y*\a k >= (Yf } ) 2 + ^P^ti (33) 

(cl) 

where Y" fc is the classical collapsing solution. Like in the homogeneous case, 
< a k \Y^\ak > does not reach the classical singularity, < Y 2 >= 0. For a 
general quantum state (29), we have the same quantum correction, so the 
expectation value is generally everywhere regular. 

In this work we use a week condition for a "quantum avoidance" of clas- 
sical singularities: given the quantum state of the system, the quantum ge- 
ometry is singular, only if its metric expectation value is. 

In the case of spherically symmetric dust universes (or stars), the expec- 
tation value < ^Ip^l^ > is everywhere regular. Semiclassically it describes 
a dust universe (or a dust star) that collapse to a minimum radius, the Planck 
radius, and does not reach the classical singularity. 

The existence of a singularity in the space-time structure is strongly re- 
lated to the "information puzzle" during the Hawking evaporation. The stan- 
dard view is that if the information is lost, it is lost "inside" the singularity. 
It may be interesting to study the Hawking radiation with a "semiclassical 
background" [9], < *|^|* >. 
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